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HIGHLIGHTS 


•  Pd-Ir  shows  high  catalytic  activity  towards  the  NaBH4  oxidation. 

•  Low  H2  generation  (<0.1  dm3  min-1)  at  potentials  more  positive  than  the  OCP. 

•  The  mechanism  of  NaBH4  oxidation  on  Pd— Ir  alloys  was  computationally  elucidated. 

•  Experimental  and  DFT  studies  agree  on  the  H2  generation  at  the  OCP. 
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A  Pd-Ir  alloy  (1:1)  coated  on  microfibrous  carbon  (11  pm  diameter)  supported  on  a  titanium  plate  was 
evaluated  as  an  electrode  for  the  anodic  oxidation  of  borohydride.  The  hydrogen  generated,  due  to  the 
parallel  reaction  of  borohydride  hydrolysis,  was  measured  during  the  electrolysis  obtaining  less  than 
0.1  cm3  min-1  H2  between  -1  and  0  V  vs.  Hg/HgO  (-0.86  and  0.14  V  vs.  SHE),  while  the  current  densities 
for  the  oxidation  of  borohydride  were  up  to  367  mA  cm-2  in  0.5  mol  dm-3  NaBH4  +  3  mol  dm-3  NaOH. 
The  low  rate  of  hydrogen  generation  suggests  that  Pd-Ir  could  be  a  promising  catalyst  for  borohydride 
oxidation.  However,  higher  rates  of  hydrogen  were  generated  at  the  open  circuit  potential,  which  is 
inconvenient  in  the  direct  borohydride  fuel  cell.  Cyclic  voltammetry  allowed  analysis  of  the  oxidation 
peaks  due  to  the  borohydride  oxidation.  To  obtain  a  further  understanding  of  the  borohydride  oxidation 
mechanism  at  Pd— Ir  electrodes,  density  functional  theory  (DFT)  was  used  to  examine  the  reaction 
mechanism  at  Pd2— Ir4( 111 )  and  Pd2— Ir2( 111 )  surfaces.  The  competition  between  borohydride  oxidation 
and  hydrogen  evolution  on  the  Pd-Ir  alloys  is  compared  with  that  on  pure  Pd(lll),  suggesting  that  the 
presence  of  Ir  favors  borohydride  oxidation  rather  than  hydrogen  evolution. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Sodium  borohydride  has  been  proposed  as  an  alternative  fuel  for 
fuel  cells  as  it  presents  many  advantages  in  comparison  with  other 
fuels:  it  is  safer  and  easier  to  transport  than  hydrogen  gas,  the 
oxidation  products  (metaborates  and  water)  are  environmentally 
friendly,  and  the  theoretical  specific  energy  density  of  the  direct 
borohydride  fuel  cell  (DBFC)  is  larger  than  other  systems, 
9.3  kW  h  kg-1.  However,  there  are  still  major  issues  that  need  to  be 
resolved  if  the  DBFC  is  going  to  advance  in  application.  These 
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include  the  lack  of  a  selective  catalyst  with  high  activity  for  direct 
borohydride  oxidation,  a  reliable  anionic  membrane,  the  instability 
of  borohydride  and  avoidance  of  its  hydrolysis.  These  aspects  limit 
the  fuel  cell  efficiency  to  less  than  30%  1,2]  of  the  theoretical 
available  energy. 

Sodium  borohydride  is  highly  soluble  in  alkaline  solution  and  is 
oxidized  at  the  anode  through  the  following  Reaction  (1): 

BH4  aq  +  80Haq  ->  B02  aq  +  6H2 Oaq  +  8e 
E°  =  -1.24  V  vs.  SHE 

However,  a  parallel  hydrolysis  Reaction  (2)  decreases  the  overall 
number  of  electrons  released  per  borohydride  molecule  to  less  than 
8,  decreasing  the  power  density  and  producing  hydrogen  gas  in  the 
anolyte: 
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BH4  aq  +  2H2Oaq  -  B02  aq  +  4H2  (2) 

Although  the  hydrogen  produced  could  be  oxidized,  the  overall 
energy  output  will  be  lower.  Hydrogen  bubbles  can  also  block  the 
electrode  surface  and  should  be  minimized  to  avoid  safety  issues. 
The  overall  Reaction  (3)  is  the  result  of  combining  Reactions  (1 )  and 
(2): 

BH4  +  xOH"  ^B02  +  (x  -  2)H20  +  (4  -  (l/2)x)H2  +  xe~  (3) 

where  x  is  the  number  of  electrons  released  per  borohydride  ion  3  ]. 

Several  catalysts  have  been  studied  for  the  direct  oxidation  of 
borohydride  ions.  Gold  and  silver  have  been  reported  to  have  good 
selectivity,  extracting  an  equivalent  8  or  almost  8  electrons  during 
oxidation  [4-6].  However,  Chatenet  et  al.,  [7]  claimed  that  the 
hydrogen  generation  rate  during  the  borohydride  oxidation  at  a 
gold  electrode  is  not  negligible  and  that  at  high  potentials 
(>0.3-0.5  V  vs.  SHE)  the  number  of  electrons  released  per  boro¬ 
hydride  molecule  is  less  than  8.  Catalysts  such  as  Pt,  Pd  and  Ni  also 
catalyze  both  oxidation  and  hydrolysis  reactions,  releasing  between 
3  and  6  electrons  per  borohydride  converted  [2,3,8-10].  The 
highest  power  density  reported  for  a  DBFC,  1.44  W  cm-2,  was  ob¬ 
tained  using  a  25  cm2  NaBH4/H202  single  cell  with  a  Pd/C  anode 
and  an  Au/C  cathode  separated  by  a  Nation  112  membrane  11  .  The 
authors  reported  that  the  effective  number  of  electrons  released  on 
a  Pd  electrode  was  less  than  8  (between  4  and  6  depending  on  the 
NaBH4  concentration)  [3,12  .  Despite  not  providing  the  full  8  elec¬ 
trons,  Pd  catalysts  lead  to  large  current  densities  and  charge  effi¬ 
ciencies  [3,12]. 

Although  Pd-based  alloys  have  not  been  tested,  alloying  elec¬ 
trodes  such  as  Pt,  Ag  or  Au  with  other  materials  has  improved  the 
performance  of  the  DBFC.  Gyenge  et  al.,  [8]  obtained  higher  peak 
currents  and  more  negative  borohydride  oxidation  potentials  by 
using  a  Pt-Au  alloy  instead  of  pure  Pt.  Alloys  such  as  Pt-Ir  and 
Pt-Ni  were  also  tested  as  working  electrodes  for  a  DBFC  [8  .  The 
anolyte  solutions  contained  2  mol  dm-3  NaBH4  in  2  mol  dm'3 
NaOH,  a  5  mg  cm-2  Pt-Ir  and  Pt-Ni  loading,  a  Nation  117  mem¬ 
brane,  and  an  02  gas  diffusion  cathode  with  4  mg  cm-2  Pt.  The  use 
of  Pt-Ir  and  Pt-Ni  catalysts  gave  the  same  power  density  of 
53  mW  cm-2  at  75  °C,  showing  the  most  active  catalytic  activity 
among  all  the  catalyst  being  tested.  When  Pt-Ir  was  used,  the 
oxidation  peaks  were  shifted  to  more  negative  potentials  than 
those  observed  with  Pt-Ni,  suggesting  the  potentially  favorable 
kinetics  on  Pt-Ir  alloy  and  yielding  the  highest  voltammetric  BH4 
oxidation  current  densities  [8  .  Kiran  et  al.,  [13]  increased  the  po¬ 
wer  density  of  a  DBFC  by  35%  by  using  an  Ir-Rh  alloy  supported  on 
carbon  compared  to  a  carbon  supported  Rh  anode. 

It  might  be  possible  to  increase  the  activity  of  pure  Pd  by 
alloying  it  with  Ir,  as  it  was  observed  by  Gyenge  et  al.,  [8]  and  Kiran 
et  al.,  [  13  ]  for  the  addition  of  Ir  to  Pt  and  Rh  to  form  Pt-Ir  and  Rh-Ir 
alloys,  respectively.  In  this  work,  cyclic  voltammetry  was  per¬ 
formed  to  analyze  the  oxidation  peaks  and  current  densities  of 
borohydride  oxidation  on  a  Pd-Ir  alloy  anode.  Also  the  hydrogen 
generation  rate  was  measured  during  the  potentiostatic  electrolysis 
of  a  borohydride  alkaline  solution  at  constant  potential.  In  order  to 
corroborate  the  experimental  work  and  obtain  a  better  under¬ 
standing  of  the  reaction  mechanism  of  the  borohydride  ions  at  the 
Pd-Ir(lll)  surface,  DFT  calculations  have  been  used  and  the 
preferred  reaction  path  is  presented.  DFT  methods  have  previously 
been  used  to  evaluate  the  borohydride  oxidation  mechanism  on 
Au(lll),  Pt(lll)  and  Au2Cui(111)  surfaces  [14-16  .  The  analysis  of 
the  Pd-Ir(lll)  elementary  surface  reaction  energetics  through  DFT 
calculations  can  help  to  corroborate  the  selectivity  of  Pd-Ir  alloys 
towards  the  oxidation  of  borohydride.  The  linear  free  energy, 
computational  hydrogen  electrode  method  is  used  to  determine 


potential  dependent  reaction  energies  [17  .  Potential  dependent 
activation  barriers  for  elementary  reaction  steps  are  determined 
using  an  approach  recently  developed  in  the  Janik  group  [16,18]. 
The  effect  of  the  Ir  content  of  in  the  Pd-Ir  alloy  is  also  analyzed  by 
considering  reaction  energetics  on  Pd2Iri( 111 )  and  Pd2Ir2( 111 ) 
surfaces. 

2.  Experimental  details 

The  electrode  material  consisted  of  Pd-Ir  nanoparticles  depos¬ 
ited  on  carbon  microfibres  of  11  pm  diameter  and  supported  on  a 
planar  titanium  foil  substrate.  The  microfibrous  carbon  catalyst 
support  was  obtained  by  a  direct  charging  electrostatic  flocking 
technique  and  was  provided  by  Dr.  C.  Patrissi  from  the  Naval  Under 
Sea  Water  Fare  Centre  [19  .  The  electrode  was  manufactured  by 
applying  a  30-100  kV  pulse  between  a  carbon  plate,  containing  the 
0.5  mm  length  and  10  pm  diameter  carbon  microfibres  (Mitsubishi 
Chemical  K637212),  and  the  titanium  foil  substrate  plate  contained 
a  carbon  adhesive  film  19].  The  density  of  the  carbon  microfibres 
on  the  titanium  foil  was  around  12,5000  fibers  per  cm-2.  Fig.  1 
shows  SEM  images  of  the  carbon  fibers  on  Ti,  a)  before  and  b)  af¬ 
ter  the  Pd-Ir  catalyst  coating.  A  catalyst  loading  of  12.3  mg  cm-2, 
with  a  composition  of  50:50  Pd:Ir,  was  obtained  by  linearly  cycling 
the  electrode  with  the  carbon  microfibres  in  a  solution  containing 

2  x  10-3  mol  dm  3  PdCl2  and  2  x  10~3  mol  dm-3  Na2IrCleH20  in 
0.2  mol  dm-3  KC1  and  0.1  mol  dm-3  HC1  at  70  °C.  Cyclic  voltam¬ 
metry  was  carried  out  at  1  mV  s_1  between  -0.15  V  and  -0.3  V  vs. 
Ag/AgCl  [19]. 

Potentiostatic  electrolysis  at  a  constant  potential  were  carried 
out  between  -0.8  and  0.4  V  vs.  Hg/HgO  (-0.66  and  0.54  V  vs.  SHE) 
in  order  to  measure  the  hydrogen  generation  rate  at  the  micro- 
fibrous  Pd-Ir  electrode.  The  accumulative  hydrogen  gas  generated 
from  the  borohydride  hydrolysis  and  the  current  densities  during 
the  electrolysis  were  monitored.  A  computer  controlled  potentio- 
stat/galvanostat  PGSTAT30  from  Autolab  (EcoChemie,  Netherlands) 
fitted  with  the  General  Purpose  Electrochemical  System  GPES  4.9 
software  was  used  to  apply  the  constant  potential.  A  three  electrode 
glass  cell  was  used  during  the  electrolysis  experiments  with  the 
working  (3  cm2  geometric  area)  and  counter  electrode  compart¬ 
ments  separated  with  a  Nation  117  ion  exchange  membrane.  The 
counter  electrode  compartment  was  filled  with  a  solution  of 

3  mol  dm-3  NaOH  while  the  working  electrode  compartment  was 
gas  tight  and  connected  to  an  inverted  50  cm3  burette  in  order  to 
measure  the  amount  of  hydrogen  gas  generated  during  borohy¬ 
dride  oxidation.  The  working  electrode  compartment  was  filled 
with  solutions  containing  3  mol  dm-3  NaOH  and  two  NaBH4  con¬ 
centrations:  0.5  mol  dm-3  or  1  mol  dm'3.  The  cyclic  voltammetry 
experiments  were  carried  out  at  10  mV  s-1  at  lower  concentrations 
of  borohydride  0.02  mol  dm-3  in  3  mol  dm'3  NaOH  in  a  cell  with  a 
smaller  Pd-Ir  electrode  of  0.283  cm2  (0.6  cm  diameter)  geometrical 


Fig.  1.  SEM  image  of  a)  carbon  fibers  deposited  by  direct  charging  electrostatic 
flocking,  b)  Pd-Ir  coated  carbon  microfibre  array  [19]. 
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area.  Both  cells  used  a  Pt  mesh  (2  cm2)  as  counter  electrode, 
separated  by  a  Nation  117  proton  membrane,  and  a  Hg/HgO  refer¬ 
ence  electrode. 

3.  Results  and  discussion 

Fig.  2  shows  the  hydrogen  generation  rate  vs.  the  electrode 
potential,  used  during  various  electrolysis  experiments.  At  each 
potential,  the  electrolysis  was  carried  out  for  60  min.  The  graph 
shows  the  curves  for  electrolysis  of  two  solutions:  a)  0.5  mol  dm  3 
and  b)  1  mol  dirT3  of  NaBH4  both  in  3  mol  dm-3  NaOH.  The  graph 
also  shows  the  temperature  measured  at  the  end  of  the  each 
experiment.  At  the  open  circuit  potential,  -0.99  V  vs.  Hg/HgO 
(-0.85  V  vs.  SHE),  i.e.  at  zero  current  density,  the  hydrogen  gen¬ 
eration  rates  were  1.3  cm3  min-1  and  0.65  cm3  min-1  for  borohy- 
dride  concentrations  of  0.5  mol  dm"3  and  1  mol  dm"3,  respectively. 
The  hydrolysis  reaction  occurs  upon  the  contact  between  the  so¬ 
dium  borohydride  solution  and  the  Pd-Ir  electrode  surface. 

At  potentials  between  -0.8  V  and  -0.4  V  vs.  Hg/HgO  (-0.66  V 
and  -0.24  V  vs.  SHE)  and  a  concentration  of  borohydride  of 
0.5  mol  dm"3  NaBH4,  the  hydrogen  generation  rate  diminished 
almost  completely  while  current  densities  from  107  mA  cm-2  to 
210  mA  cm-2  were  obtained.  At  more  positive  potentials  than  -0.2  V 
vs.  Hg/HgO  (-0.06  V  vs.  SHE),  the  gas  generation  rate  slightly  started 
to  increase,  reaching  0.96  cm3  min-1,  at  0  V  vs.  Hg/HgO  (0.14  V  vs. 
SHE).  In  contrast,  at  1  mol  dm  3  NaBH4  a  lower  hydrogen  generation 
rate  was  measured  at  the  open  circuit  potential  and  higher  hydrogen 
generation  rate  was  obtained  at  all  the  other  measured  potentials. 
This  could  be  related  to  the  Pd  oxidation  states  [3];  at  low  concen¬ 
trations  of  borohydride  ions  the  electrode  surface  remains  in  a 
passive  state,  the  Pd  on  the  surface  oxidizes  and  the  borohydride 
oxidation  decreases.  Lower  current  densities  and  lower  hydrogen 
generation  are  expected  for  lower  concentrations  of  borohydride.  At 
high  concentrations  of  borohydride  ions,  the  surface  oxide  or  hy¬ 
droxide  might  be  reduced  by  the  borohydride  and  become  active 
again  towards  the  oxidation  of  borohydride  ions  [3],  also  leading  to 
more  hydrogen  evolution.  At  -0.4  V  vs.  Hg/HgO  (-0.24  V  vs.  SHE), 
the  hydrogen  generation  rate  increased  from  0.65  cm3  min-1  to 
1.28  cm3  min-1,  when  the  concentration  of  borohydride  was 
increased  from  0.5  mol  dm"3  to  1  mol  dm-3,  respectively.  This 
agrees  with  the  suggestion  of  Liu  et  al.,  [3]  that  at  palladium  elec¬ 
trodes,  higher  efficiency  (6e-)  is  obtained  using  relatively  low  con¬ 
centrations  of  borohydride,  whereas  for  concentrations  higher  than 
1  mol  dm-3  the  efficiency  decreases  (4e-). 


Fig.  2.  Gas  generation  rate  vs.  potential  applied  at  the  Pd/Ir  coated  microfibrous  carbon 
electrode  (3  cm2),  a)  A  1  mol  dm'3  NaBH4  in  3  mol  dm-3  NaOH,  b)  ▲  0.5  mol  dm-3 
NaBH4  in  3  mol  dm  3  NaOH. 


During  electrolysis  of  the  electrolyte  containing  1  mol  dm"3 
NaBH4,  the  temperature  increased  from  23  °C  at  the  beginning  of 
the  experiment  to  50,  70  and  78  °C  at  -0.26,  -0.06  and  0.14  V  vs. 
SHE,  respectively.  The  increase  in  temperature  is  due  to  the 
exothermic  oxidation  and  hydrolysis  reactions  that  clearly 
increasing  the  hydrogen  generation  rate  [20,21]  as  can  be  seen  in 
curve  a)  of  Fig.  2.  At  a  0.5  mol  dm-3  NaBH4  concentration,  the 
temperature  increased  from  23  °C  to  30  °C  at  all  the  applied  po¬ 
tentials.  The  lower  temperature  increase  during  the  experiments 
carried  out  with  0.5  mol  dm-3  NaBH4  was  due  to  the  near  absence 
of  exothermic  hydrolysis  compared  to  that  occurring  with 
1  mol  dm-3  NaBH4. 

Fig.  3  shows  a  comparison  of  the  hydrogen  generation  rate  vs. 
current  for  1  mol  dm-3  and  0.5  mol  dm-3  NaBH4  in  3  mol  dm"3 
NaOH,  curves  a)  and  b),  respectively  during  constant  potential 
electrolysis  at  potential  values  between  -0.8  V  and  0  V  vs.  Hg/HgO 
(-0.66  V  and  0.14  V  vs.  SHE).  The  current  density  plotted  in  Fig.  3 
was  measured  during  the  electrolysis  at  constant  potential  and 
the  hydrogen  evolution  rate  at  each  point  corresponds  to  the  po¬ 
tentials  of  the  curves  a)  and  b)  in  Fig.  2.  The  figure  also  compares 
the  results  obtained  in  this  work  with  curves  c)  and  d),  reported  by 
Liu  et  al.,  [3]  using  a  Pd/C  electrode  having  an  8  mg  cm-2  catalyst 
loading  and  solutions  of  0.52  mol  dm"3  in  2  mol  dm-3  NaOH  and 
1  mol  dm"3  NaBH4  in  2  mol  dm"3  NaOH,  respectively.  All  the  curves 
show  similar  trend:  hydrogen  is  generated  at  zero  current,  followed 
by  a  decrease  in  the  hydrogen  generation  rate  between 
100  mA  cm-2  and  400  mA  cm-2  and  an  increase  of  the  hydrogen 
generation  rate  at  current  densities  from  400  mA  cm-2  to  higher 
values.  Liu  et  al.,  obtained  a  hydrogen  generation  rate  of 
4  cm3  min-1  at  100  mA  cm-2  from  the  solution  containing 
1  mol  dm-3  NaBH4  in  2  mol  dm"3  NaOH  at  a  potential  of  -0.5  V  vs. 
SHE  (curve  d).  At  the  same  current  density,  a  much  lower  hydrogen 
generation  rate,  0.1  cm3  min-1,  was  measured  in  this  work  using  a 
Pd-Ir  alloy  and  the  same  concentration  of  borohydride  in 
3  mol  dm-3  NaOH.  At  300  mA  cm-2  Liu  et  al.,  obtained  a  hydrogen 
generation  rate  of  nearly  2.28  cm3  min-1  and  1.57  cm3  min-1,  for 
1  mol  dm-3  and  0.5  mol  dm"3  NaBH4  in  2  mol  dm-3  NaOH, 
respectively  whereas  values  below  0.5  cm3  min-1  were  obtained 
using  the  Pd/Ir  coated  microfibrous  carbon  electrode  at  the  same 
current  density  in  this  work.  The  large  difference  in  the  hydrogen 
generation  rate  obtained  by  Liu  et  al.,  [3]  compared  to  that  obtained 


Fig.  3.  Gas  generation  rate  vs.  current  from  electrolysis  at  constant  potential:  a)  ▲  Pd/ 
Ir  coated  microfibrous  carbon  electrode  using  a  solution  of  0.5  mol  dm  3  NaBH4  in 
3  mol  dm-3  NaOH,  b)  A  Pd/Ir  coated  microfibrous  carbon  electrode  and  1  mol  dm'3 
NaBH4  in  3  mol  dm'3  NaOH,  c)~10  wt.  %  Pd/C  and  0.52  mol  dm  3  NaBH4  in  2  mol  dm'3 
NaOH  [3  ,  d)  •  10  wt.  %  Pd/C  and  1.02  mol  dm'3  NaBH4  in  2  mol  dm'3  NaOH  [3]. 
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in  this  work  is  attributed  to  the  difference  in  concentration  of 
NaOH,  the  presence  of  Ir  and  a  higher  loading  of  catalyst,  which 
improved  the  activity  by  increasing  the  oxidation  of  borohydride 
and  decreasing  the  rate  of  hydrogen  generation. 

Fig.  4  shows  the  apparent  number  of  electrons  released  per 
borohydride  molecule  converted  during  the  electrolysis  at  constant 
potential.  This  value  can  be  calculated  according  to  Equation  (4), 
proposed  by  Wang  et  al.,  [22]. 


where  /  is  the  measured  current  and  Iq  is  the  current  due  to 
oxidation  hydrogen  generated  during  the  hydrolysis  of  borohydride 
ions  (calculated  using  the  Faraday's  law  and  considering  a  2e~ 
transfer).  For  a  complete  borohydride  oxidation,  napp  =  8,  there  is 
no  hydrogen  generation  and  Ic  will  be  zero.  According  to  Fig.  4,  the 
number  of  electrons  transferred  during  the  electrolysis  of 
0.5  mol  dm-3  NaBFLi  in  3  mol  dm-3  NaOH  was  almost  8  for  po¬ 
tentials  between  -0.6  and  -0.4  V  vs.  Hg/HgO  (-0.46  and  -0.26  V 
vs.  SHE),  whereas  at  1  mol  dm-3  NaBH4  in  3  mol  dm-3  NaOH  the 
number  of  electrons  decreased  towards  more  positive  potentials 
from  6e_  at  -0.8  V  vs.  Hg/HgO  (-0.66  V  vs.  SHE)  to  less  than  2e~ 
at  -0.2  V  vs.  Hg/HgO  (-0.06  V  vs.  SHE). 

The  cyclic  voltammogram  shown  in  Fig.  5  is  for  a  0.02  mol  dm  3 
BH4  ions  concentration  at  the  Pd-Ir  coated  microfibrous  carbon 
electrode  (3  cm2)  in  a  three-electrode  electrochemical  cell.  The 
Figure  also  shows  the  CV  of  a  blank  solution  (3  mol  dm  3  NaOH  in 
the  absence  of  borohydride)  for  comparison.  A  potential  peak,  /, 
obtained  at  -0.95  V  vs.  Hg/HgO  (-0.81  V  vs.  SHE)  is  attributed  to 
the  oxidation  of  the  hydrogen  adsorbed  on  the  electrode  surface 
[23].  In  the  presence  of  the  borohydride  ions,  the  peak  is  noticeably 
increased  due  to  the  presence  of  hydrogen  formed  during  the  hy¬ 
drolysis  of  borohydride  at  the  open  circuit  potential  (-0.99  V  vs. 
Hg/HgO  or  -0.85  V  vs.  SHE).  The  second  peak  II,  also  observed  in  the 
blank  solution,  at  -0.71  V  vs.  Hg/HgO  (-0.57  V  vs.  SHE),  is  sur¬ 
rounded  by  small  shoulders  indicating  a  complex  oxidation  process 
of  palladium  24].  This  peak  is  attributed  to  the  adsorption  of  OH~ 
in  the  initial  stage  of  the  palladium  oxide  formation  and  partially 
overlaps  the  hydrogen  oxidation  peak  [23  .  The  small  rise  III,  which 
emerges  from  peak  II,  might  be  due  to  the  second  stage  of  forma¬ 
tion  of  the  palladium  (II)  oxide  layer  on  the  surface  of  the  catalyst 


Fig.  4.  Apparent  number  of  electrons  released  during  oxidation  vs.  electrode  potential 
from  solutions  containing:  a)  0.5  mol  dm-3  NaBH4  in  3  mol  dm  3  NaOH,  b)  1  mol  dm-3 
NaBH4  in  3  mol  dm  3  NaOH  at  23  °C  using  an  Pd-Ir  electrode. 


Fig.  5.  Cyclic  voltammogram  at  10  mV  s^1  at  the  Pd/Ir  coated  microfibrous  carbon 
electrode  (3  cm2).  Solutions  containing  0.02  mol  dm  3  NaBH4  in  3  mol  drrr3  NaOH 
(black  line)  and  3  mol  dm  3  NaOH  (broken  line),  were  used. 

[25].  In  the  presence  of  borohydride,  the  increase  in  current  den¬ 
sities  of  peaks  II  and  III  is  probably  due  to  a  mixed  oxidation  of 
palladium,  borohydride  ions  and  possibly  intermediate  species 
formed  during  borohydride  oxidation.  The  reduction  peak  IV 
at  -0.79  V  vs.  Hg/HgO  (-0.65  V  vs.  SHE)  on  the  reverse  scan, 
observed  both  in  the  presence  of  borohydride  and  in  the  blank 
solutions  can  be  attributed  to  the  reduction  of  Pd(II)  oxide  to  Pd(0). 
On  the  reverse  scan  of  the  borohydride  solutions  CV  an  additional 
oxidation  peal  V  is  observed  on  the  reverse  scan  at  -0.17  V  vs.  Hg / 
HgO  (-0.03  V  vs.  SHE),  which  can  be  attributed  to  an  increased  rate 
of  borohydride  oxidation  on  the  Ir  atoms  or  on  the  more  active 
metal  formed  during  the  reduction  of  PdO.  Kiran  et  al.,  [13]  ob¬ 
tained  an  oxidation  peak  at  -0.43  V  vs.  Hg/HgO  (-0.29  V  vs.  SHE) 
when  they  carried  out  the  CV  of  the  0.1  mol  dm-3  NaBH4  in 
1  mol  dm-3  NaOH  using  an  Ir/C  anode. 

The  results  shown  above  suggested  that  the  Pd-Ir  catalyst  has 
activity  towards  borohydride  oxidation,  resulting  in  current  den¬ 
sities  up  to  367  mA  cm-2  during  the  electrolysis  at  concentrations 
of  0.5  mol  dm-3  of  BH4.  Meanwhile,  the  hydrogen  generation  rate 
was  negligible  from  -0.8  to  0.2  V  vs.  Hg/HgO  (-0.66  to  -0.06  V  vs. 
SHE).  At  higher  BH4  concentrations,  the  hydrogen  generation 
increased  from  -0.4  V  vs.  Hg/HgO  (-0.26  V  vs.  SHE)  together  with 
the  temperature  of  the  reaction.  However,  the  hydrogen  generated 
at  the  OCP  was  significant  for  both  concentrations  and  might  be  a 
problem  for  the  cell  performance.  The  results  also  showed  that  the 
presence  of  Ir  improved  the  catalytic  activity  of  Pd.  In  order  to 
corroborate  the  positive  effect  of  Ir  in  the  Pd-Ir  alloy  and  to  obtain  a 
better  understanding  of  the  reaction  mechanism  of  borohydride 
oxidation  on  Pd-Ir  surfaces,  DFT  calculations  of  this  system  are 
presented  in  the  following  section. 

4.  Computational  methods 

DFT  calculations  were  performed  to  elucidate  the  reaction 
mechanism  of  borohydride  oxidation  on  Pd-Ir(lll)  surfaces.  The 
effect  of  the  Ir  concentration  in  the  Pd-Ir  alloy  was  analyzed  by 
studying  the  oxidation  mechanism  on  two  different  bimetallic 
Pd2lri( 111 )  and  Pd2lr2(lll),  surfaces.  The  optimum  reaction  path 
was  obtained  by  investigating  the  elementary  surface  reaction 
energies  and  postulating  the  lowest  energy  path  taking  into  ac¬ 
count  the  activation  energy  for  each  electron  transferred.  As  a  first 
approximation,  the  energy  of  a  transferred  electron  and  the  ion 
adsorption  energies  were  made  linearly  dependent  on  the 
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electrochemical  potential  and  the  energy  of  the  adsorbed  species 
was  calculated  as  their  free  energy  14,17,26  . 

All  calculations  were  performed  using  the  ab  initio  total  energy 
and  molecular-dynamics  Vienna  ab  initio  simulation  program 
(VASP)  developed  at  the  Institute  for  Material  Physics  at  the  Uni¬ 
versity  of  Vienna  [27—29].  The  projected  augmented  wave  method 
[30,31  and  a  plane  wave  basis  set  were  employed  using  an  energy 
cut  off  of  450  eV.  The  generalized  gradient  approximation  (GGA) 
based  on  that  of  Perdew-Wang  32]  was  used  to  calculate  the  ex¬ 
change  and  correlation  energies.  To  optimize  the  structure,  a 
3x3x1  Monkhorst-Pack  grid  [33]  was  used,  followed  by  a 
5x5x1  grid  single-point  calculation  to  give  the  total  energy. 

The  Pd2lri  and  Pd2lr2  bimetallic  surfaces  were  modeled  using  a 
3  layer  slab  and  a  3  x  3  surface  cell  and  a  3  layer  slab  with  a  4  x  4 
surface  cell,  respectively;  with  the  two  bottom  layers  constrained  to 
fee  lattice  positions  and  the  two  top  layers  relaxed.  A  vacuum  of 
13  A  was  inserted  between  the  periodic  slabs  as  a  representation  of 
the  electrode  surface.  The  borohydride  molecule  and  all  possible 
borohydride  oxidation  intermediates  were  optimized  on  both 
Pd2lri  and  Pd2lr2  surfaces. 

The  adsorption  energies  of  BH4  and  borohydride  oxidation  in¬ 
termediates  were  calculated  at  the  optimal  adsorption  configura¬ 
tions  and  were  corrected  for  spurious  slab-to-slab  dipole 
interactions  along  the  surface  normal  direction.  The  reaction  en¬ 
ergy  and  activation  energy  for  each  step  of  the  oxidation  of  boro¬ 
hydride  were  calculated  using  the  lowest  energy  adsorbed  states 
[15,16  .  The  harmonic  vibrational  modes  were  calculated  to  deter¬ 
mine  the  zero  point  vibrational  (ZPVE)  corrections  to  the  total  en¬ 
ergy  of  adsorbates.  Anderson  and  Kang  [26]  and  Norskov  et  al.,  [17] 
have  used  similar  methods,  with  various  perturbations,  to  study  a 
number  of  electrocatalytic  reactions.  The  free  energies  of  aqueous 
phase  molecules  and  ions  such  as  BH4,  BO2 ,  H2O  and  OH-  were 
previously  reported  by  Rostamikia  et  al.,  [14],  who  calculated  those 
values  by  adding  ZPVE,  pressure  volume  (PV)  and  entropic  terms  to 
the  isolated  gas  phase  molecule  energy  (EDFr).  Experimental  sol¬ 
vation  free  energies  and  concentration  adjustments  were  included 
to  consider  solution  conditions,  assuming  an  ideal  solution. 

The  activation  energy  for  each  step  of  the  borohydride  oxidation 
reaction  on  Pd-Ir  surfaces  was  calculated  by  locating  transition  states 
between  the  lowest  energy  adsorbed  states,  presuming  that  each  step 
follows  the  minimum  energy  path  [15,16  .  The  transition  states  of 
elementary  reactions  on  the  surfaces  were  isolated  with  the  climbing 
image  nudged  elastic  band  method  (CI-NEB)  [34  ,  which  finds  saddle 
points  and  minimum  energy  paths  between  the  known  reactants  and 
products.  Elementary  reactions  for  which  barriers  were  computed 
include  B-H  dissociation  reaction  and  B-OH  formation  reactions.  A 
number  of  equally  spaced  images  (usually  4  images  and  occasionally 


8)  between  the  initial  and  final  states  were  optimized,  with  reaction- 
tangent  forces  and  atomic  forces  meeting  optimization  criteria  and  a 
single  harmonic  imaginary  frequency  along  the  reaction  coordinate. 
The  highest  energy  image  corresponds  to  the  transition  state,  which 
gives  the  activation  energy  of  the  reaction. 

The  B-H  dissociation  or  B-OH  formation  activation  barriers 
were  extrapolated,  using  a  Butler-Volmer  formalism  and  an 
approach  recently  developed  in  the  Janik  group  16],  in  order  to 
make  them  potential  dependent.  A  symmetry  factor  is  used  to 
linearly  relate  the  activation  barriers  to  the  equivalent  electro¬ 
chemical  reaction  15]: 

K(E)  =  ^-P(E-E0)  (5) 

where  0^  is  the  activation  barrier  for  the  chemical  B-H  dissocia¬ 
tion  or  B-OH  formation  step,  and  (3  is  the  symmetry  factor  which 
will  be  considered  0.5,  assuming  a  symmetric  reaction  and  in  the 
range  of  most  of  the  electrochemical  reactions  [35  .  E0  is  the 
equilibrium  potential  that  converts  the  ion  involved  in  the  elec¬ 
trochemical  step  into  the  equivalent  surface  species  used  in 
locating  the  transition  state  (H+  or  OH~  converting  to  H  or  OH" 
with  electron  transfer). 

5.  DFT  calculation  results  and  discussions 

5.1.  Borohydride  ion  adsorption  over  Pd2—lr\( W)  and  Pd2-Ir2( 111) 

Fig.  6(a)  and  (b)  shows  the  optimum  configuration  for  the  mo¬ 
lecular  adsorption  of  borohydride  ions  on  Pd2lri( 111 )  and 
Pd2lr2(lll)  surfaces,  respectively.  The  BH4  species  adsorbs  on  the 
Pd2 — Iri(lll )  slab  with  the  boron  atom  above  the  Pd-Pd  bridge  and 
two  H  atoms  on  Pd  atop  sites  as  shown  in  Fig.  6(a).  This  was  the  only 
stable  configuration  found  for  a  molecularly  adsorbed  (BH4)  on 
Pd2 — Iri(lll ).  The  activation  energy  for  the  decomposition  of  BH4  to 
BH3  +  H*  is  negligible,  dissociative  adsorption  occurs  spontane¬ 
ously  and  is  thermodynamically  favorable  (AGads  <  0),  producing 
BH3  +  H*.  The  adsorption  reaction  can  be  written  as: 

BH4  +  *  — >  (BH3  +  H*)  +  e~  (6) 

where  *  represents  a  bare  surface.  As  shown  in  Fig.  6(b),  the 
borohydride  ion  is  adsorbed  on  the  Pd2lr2(lll)  surface  with  the 
boron  atom  above  the  hollow  site,  two  hydrogen  atoms  on  the  Pd 
atoms  atop  sites  and  one  hydrogen  atom  on  the  Ir  atom  atop  site. 
The  adsorption  reaction  of  the  borohydride  ion  on  the  Pd2lr2(lll) 
surface  corresponds  to  Reaction  (7): 


Fig.  6.  BH4  molecule  adsorbed  on:  a)  Pd2-Iri(lll)  with  two  H  atoms  on  Pd  atoms  and  the  other  two  H  atoms  on  atop  sites,  b)  Pd2 — Ir2(lll ).  The  dark  blue  atoms  correspond  to  Ir, 
the  light  blue  slightly  larger  atoms  correspond  to  Pd,  the  white  atoms  are  H  atoms  and  the  pink  atoms  are  B.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.) 
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BH4(aq)  +  *-BH;  +  e-  (7) 

The  free  energy  of  adsorption  of  borohydride  ions  on  the 
Pd2 — Iri( 111 )  and  Pd2-Ir2(lll)  surfaces  can  be  calculated  using 
Equations  (8)  and  (9): 


^^ads 


cbh;+h*  +  eE  -  C*  -  GBH4aq 

(8) 

gbh;  +eE-G*  -  CBH4aq 

(9) 

where  GBH4a  is  the  free  energy  of  the  borohydride  ions  in  aqueous 
solution  15],  G*  is  the  free  energy  of  the  bare  surface,  GBH*+H*  is  the 
free  energy  of  the  dissociatively  adsorbed  borohydride  species, 
Gbh*  is  the  free  energy  of  the  adsorbed  molecule,  e~  is  the  charge  of 
an  electron  and  E  is  the  absolute  (vacuum  reference)  electrode 
potential.  The  absolute  potential,  E ,  can  be  referred  to  a  standard 
hydrogen  electrode  using  the  equation: 


^SHE  =  E  ~  £ref,SHE  =  £  ~  4.6  (10) 

where  Eshe  is  the  potential  vs.  the  standard  hydrogen  electrode 
(SHE).  Note  that  Eshe  is  what  we  referred  to  as  E  vs.  SHE  in  the 
experimental  sections  and  in  Equation  (5).  The  SHE  potential  has  an 
average  value  of  4.6  V  referenced  to  the  absolute  vacuum  [36  ,  and 
the  value  of  4.6  V  is  consistent  with  that  calculated  by  Rostamikia 
et  al.,  15]  using  DFT  methods. 

Fig.  7  plots  the  adsorption  energy  of  borohydride  ions  over 
Pd(lll),  Pd2— Iri( 111 )  and  Pd2-Ir2(lll)  vs.  potential.  One  adsorbate 
molecule  per  nine  surface  atoms  was  considered  on  the  Pd(lll )  and 
on  the  Pd2lri  surfaces  and  one  adsorbate  molecule  per  sixteen 
surface  atoms  on  the  Pd2lr2  surface.  The  adsorption  free  energy  of 
BH4  over  the  Pd(lll),  Pd2lr4( 111 )  and  Pd2lr2(lll)  surfaces 
were  -1.77  eV,  -1.12  eV  and  -0.89  eV,  respectively  at  -0.5  V  vs. 
SHE.  The  increased  fraction  of  Ir  in  the  bimetallic  increases  the 
adsorption  energy,  making  it  less  exergonic  (releases  less  energy). 

The  adsorption  of  borohydride  over  the  Pd-Ir(lll)  surfaces  is 
stronger  than  that  on  Au(lll)  and  weaker  than  on  Pt(lll).  The 
adsorption  of  borohydride  ions  at  Au(lll)  starts  being  favorable 
(AG  <  0)  above  0.05  V  vs.  SHE  15],  whereas  at  Pd-Ir(lll)  the 
adsorption  of  borohydride  is  favorable  at  any  potential  above  -1  V  vs. 
SHE.  Pd — Ir(lll )  is  therefore,  active  towards  this  initial  oxidation  step 
at  more  negative  potentials  than  Au(  111 ).  The  borohydride  adsorption 
energy  on  Pt(lll)  was  -1.85  eV,  at  -0.5  V  vs.  SHE,  which  causes 
dissociation  to  BH*  +  3H*  [14].  At  Pd2lri(lll),  the  borohydride 


adsorption  is  less  negative  than  on  Pt(  111 )  and  the  molecule  is  spon¬ 
taneously  dehydrated  to  BH3  +  H*,  whereas  at  Pd2lr2  the  borohydride 
molecule  was  stable  and  optimized  as  BH4.  The  presence  of  Ir  makes 
the  borohydride  adsorption  energy  more  positive  and  less  likely  to  be 
dissociative  [14,37].  Rostamikia  etal.,  37  previously  suggested  that  a 
BH4  adsorption  energy  intermediate  between  those  at  Au  and  Pt,  with 
less  tendency  to  dissociate,  would  allow  for  a  better  balance  of 
increased  oxidation  activity  (compared  to  Au)  and  increased  selec¬ 
tivity  to  direct  oxidation  over  hydrolysis  (compared  to  Pt). 

5.2.  Elementary  BH4  oxidation  energiesat  Pd2lr\(lll )  and 
Pd2lr2( III)  surfaces 

The  free  energies  of  the  electro-oxidation  steps  that  involve  the 
surface  species  can  be  calculated  using  the  following  equation: 

R*  +  xOH-q  ->  O*  +  yH20  +  ne~  (11) 

where  R  is  the  reactant  adsorbed  on  the  surface,  O  is  the  adsorbed 
oxidized  product  and  n  is  the  number  of  electrons  exchanged.  The 
free  energy  of  Reaction  (11)  can  be  calculated  as: 

AGrXn  =  Gp*  +yGH20aq  —  neU  —  Gr*  —  xGqh“  (12) 

The  reaction  energy  diagrams  for  the  borohydride  oxidation 
over  Pd2 — Ir4( 111 )  and  Pd2-Ir2(lll)  surfaces  were  evaluated 
considering  the  relative  free  energies  of  all  adsorbed  species  with 
respect  to  the  free  energy  of  the  initial  aqueous  borohydride  ion. 
The  relative  energy  of  each  species  was  calculated  by  combining 
Equations  (7)  or  (8)  and  (12),  leading  to  Equation  (13)  [14]. 

Relative  energy  of  intermediate  species 

BOxH*  =  G(BOxH*)  +  (4  -  2x  +  y)C0H-  -  (4  +  2 x-y)eU  (13) 

+  (4  +  x  -  j/)Gh2 oaq  -  CBH-  -  8G0Hiq 

where  0  <  x  <  2  and  0  <  y  <  4.  The  potential  dependent  activation 
barriers  were  also  calculated  for  each  electro-oxidation  step.  As  an 
example,  Fig.  8  shows  the  initial  state  (a),  transition  state  (b)  and 
final  state  (c)  of  the  dehydrogenation  reaction  from  BH3  to  BH2  +  H* 
on  Pd2Iri.  The  activation  barrier  for  this  reaction  was  0.09  eV  (kj 
mol-1)  at  -0.5  V  vs.  SHE.  The  dissociating  B-H  bond  stretches  from 
1.2  A  in  the  initial  state  to  1.35  A  at  the  transition  state.  The  initial, 
transition  and  final  states  of  the  same  reaction  on  Pd2Ir2  are  illus¬ 
trated  in  Fig.  9(a)— (c),  respectively.  In  this  case,  at  -0.5  V  vs.  SHE, 
the  activation  barrier  was  0.36  eV  (kj  mol-1)  and  the  transition 
state  much  closer  to  the  final  state,  which  indicates  that  this  re¬ 
action  step  is  more  difficult  to  perform  on  Pd2Ir2  than  on  Pd2Iri. 

The  activation  energies  for  each  electro-oxidation  step  and 
relative  energy  of  each  species  are  used  to  evaluate  the  most  likely 
oxidation  reaction  mechanism.  Fig.  10  shows  the  relative  energies 
of  all  possible  borohydride  oxidation  intermediates  on  Pd2Iri( 111 ) 
at  -0.5  V  vs.  SHE.  In  this  diagram,  each  step  to  the  right  is  a  le- 
oxidation  reaction.  Steps  involving  B-H  dissociation  include  the 
conversion  of  a  hydroxyl  species  to  a  water  molecule,  for  example 
the  BH3  to  BH2O*  step  is,  in  complete  notation: 

BH;  +  OH-q-BH;  +  H2Oaq  +  e~  (14) 

For  reactions  including  B-OH  formation,  the  complete  reaction 
is,  for  example: 

BOH*  +OH“q— >B(OH)2  +  e~  (15) 


Fig.  7.  dft  calculated  adsorption  free  energy  of  BH4aq  over  Pd(in),  Pd2iri( ill )  and  The  most  favorable  mechanism  is  highlighted  in  the  diagram 

Pd2ir2(  111 )  using  vacuum  slab  model  (7=  298  k,  [BH^aq]  =  0.03  mol  dm  3).  with  a  red  straight  line.  The  broken  and  dotted  lines  represents 
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Fig.  8.  BH3  dehydrogenation  to  BH^+Hf  over  Pd2Ir1(lll)  surface:  a)  initial  state,  b)  transition  state,  c)  final  state.  The  dark  blue  atoms  correspond  to  Ir,  the  light  blue  slightly  larger 
atoms  correspond  to  Pd,  the  white  atoms  are  H  atoms  and  the  pink  atoms  are  B.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web 
version  of  this  article.) 


Fig.  9.  BH3  dehydrogenation  to  BH2+H*  over  Pd2lr2(lll)  surface:  a)  initial  state,  b)  transition  state,  c)  final  state.  The  dark  blue  atoms  correspond  to  Ir,  the  light  blue  big  atoms 
correspond  to  Pd,  the  pink  atoms  are  B,  the  red  are  0  atoms  and  the  white  are  H  atoms.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the 
web  version  of  this  article.) 


the  other  possible  mechanisms  calculated  with  DFT,  the  most 
favorable  mechanism  was  determined  assuming  that  the  reaction 
follows  the  lowest  activation  energies  and  the  intermediate 
molecules  with  the  lowest  relative  energies  are  involved  in  the 
reaction.  The  values  of  the  activation  energy  (eV)  for  each  possible 
step  are  shown  in  Fig.  10,  between  the  species.  The  8e~  transfer  is 
favorable  since  the  reaction  energy  is  downhill  and  each  inter¬ 
mediate  species  has  a  lower  relative  energy  than  the  previous 
one.  In  general,  B-H  breaking  elementary  steps  are  more 


favorable  than  B-OH  bond  forming  steps.  B-H  dissociation  pro¬ 
ceeds  until  the  formation  of  BH*.  The  first  reaction  step,  from  BH4 
to  BH3,  is  barrier-less  and  the  energy  of  BH3  is  lower  than  that  of 
BH4.  The  second  step  considers  two  possibilities:  1)  BH3  might 
react  with  OH"  to  form  BH3OH*  or  2)  it  might  dehydrogenate  to 
form  BH2.  The  latter  is  more  favorable  since  the  reaction  product 
BH2  has  a  lower  relative  energy  and  the  activation  barrier  for  BH2 
formation  is  0.09  eV,  compared  with  0.88  eV  for  BH3OH* 
formation. 


Fig.  10.  Reaction  free  energy  diagram  for  borohydride  oxidation  over  Pd2Iri( 111 )  at  -0.5  V  vs.  SHE,  [BH4  ]  =  0.02  mol  dm-3,  [OH-]  =  2  mol  dm  3,  T  =  298  K.  The  red  line  represents 
the  most  favorable  mechanism  of  borohydride  oxidation  according  to  the  relative  energies  and  activation  barriers  calculated  from  DFT.  The  broken  and  dashed  lines  are  other 
possible  mechanisms  of  borohydride  oxidation.  Numbers  on  the  plot  label  the  activation  barriers  in  eV.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  this  article.) 
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Following  B-H  dissociation  steps  to  BH",  reaction  with  OH" 
forms  BHOH*,  followed  by  a  low  energy  barrier  process  forming 
BOH".  From  that  point  forward,  the  activation  barriers  are  high 
enough  to  conclude  that  B-OH  bonds  do  not  form  at  a  substantial 
rate  at  room  temperature  and  at  -0.5  V  vs.  SHE.  Several  transition 
state  searches,  considering  different  reaction  trajectories  were 
performed  to  search  for  paths  with  lower  activation  barriers, 
allowing  reaching  a  final  8e-  product.  However,  the  lowest  acti¬ 
vation  energies  found  are  reported  herein  and  suggest  that  B-OH 
formation  will  be  challenging  under  these  conditions,  on  the  single 
crystal  surfaces. 

There  are  two  possible  interpretations  of  these  results: 

1 )  The  reactions  of  intermediate  products  with  hydroxyl  ions  limit 
the  overall  oxidation  reaction.  The  barriers  from  BOH  to  further 
oxidation  will  limit  the  overall  oxidation  rate  and  are  possibly 
too  high  for  the  reaction  to  go  forward,  leaving  BOH  as  the  final 
oxidation  product.  Being  a  stable  species,  BOH,  might  cover  the 
catalyst  surface  and  poison  it,  reducing  the  catalyst  activity  to¬ 
wards  the  borohydride  oxidation.  The  number  of  electrons 
released  per  BH4  converted  on  pure  Pd  and  pure  Ir  have  been 
reported  to  be  less  than  8,  and  thus  in  a  Pd-Ir  alloy  an  8  electron 
transfer  was  not  expected.  The  number  of  electrons  released, 
according  to  the  mechanism  of  oxidation  highlighted  in  red  in 
Fig.  10  on  Pd2lri,  would  be  5  terminating  with  BOH  as  the  final 
product. 

2)  The  barriers  for  B-OH  forming  reactions  from  the  DFT  approach 
used  are  generally  too  high  due  to  the  lack  of  surface  species 
solvation  in  the  surface  model.  The  complete  oxidation  reaction 
is  favorable  since  all  the  intermediate  products  are  lower  in 
energy  than  the  previous  one,  however,  the  high  barriers  hinder 
the  reaction  to  reach  the  final  product  to  release  8e-.  Solvation 
would  lower  the  barrier  in  three  ways:  water  molecules  could 
more  strongly  solvate  the  transition  state  than  initial  state,  OH" 
could  be  shuttled  through  a  H20  molecule  (H  from  H20  transfers 
to  OH"  while  a  B-OH  bond  is  formed  with  the  O  of  a  water 
molecule),  and/or  the  E0  for  OH"  formation  could  be  lowered  due 
to  solvation  of  OH*.  Nie  et  al.,  [18]  studied  the  kinetics  of 
elementary  steps  in  C02  electro-reduction,  finding  that  for  all 
O-H  bond  formation  reactions,  water  solvation  in  ion  shuttling 
through  H20  lowered  activation  barriers.  C-H  formation  re¬ 
actions  did  not  experience  similar  stabilization.  It  is  possible 
that  lower  activation  barriers  can  be  found  when  considering 


water  molecules  close  to  the  borohydride  oxidation  intermedi¬ 
ate  steps,  leading  to  easier  B-OH  bond  formation  and  being  able 
to  reach  the  final  oxidation  product,  hydrated  metaborate,  and 
releasing  8e~.  The  lack  of  solvation  of  OH"  certainly  leads  to  E0 
being  unreasonably  high,  which  raises  barriers  at  -0.5  V  vs.  SHE. 
Hydrogen  bonding  to  OH"  is  expected  to  lower  E0  by  approxi¬ 
mately  0.4  eV,  suggesting  B-OH  formation  barriers  should  be 
lower  than  those  reported  herein  by  at  least  0.4  eV. 

Fig.  11  shows  all  the  reaction  free  energy  diagram  for  the 
oxidation  of  borohydride  ions  on  Pd2Ir2( 111 ).  The  preferable 
mechanism  of  reaction  is  highlighted  with  red  lines.  The  activation 
energies  of  each  step  are  shown  in  the  figure  for  each  possible  step. 
The  reaction  starts  with  borohydride  adsorption,  which  was  opti¬ 
mized  as  a  stable  BH4.  The  transition  from  BH4  to  BH3  +  H*  was 
barrier-less,  meaning  that  BH3  formation  will  spontaneously  occur. 
The  next  oxidation  step  following  the  red  line,  oxidation  of  BH3  had 
only  one  possible  product,  BH2,  since  the  intermediate  molecule 
BH3OH"  was  unstable  and  dissociated  during  the  optimization 
process.  An  activation  barrier  of  0.36  eV  allows  the  BH3  molecule  to 
form  BH2,  followed  by  the  formation  of  BH*,  which  would  tend  to 
overcome  0.45  eV  to  form  a  more  stable  product  (BHOH*),  rather 
than  overcome  a  slightly  higher  barrier  0.49  eV  to  form  a  less  stable 
B",  that  follows  the  dotted  line.  From  BHOH"  a  low  barrier  of  0.19  eV 
led  to  the  formation  of  BOH".  Similar  to  the  mechanism  of  boro¬ 
hydride  oxidation  on  Pd2Iri( 111 ),  from  this  point  the  activation 
barriers  to  find  the  final  product  (metaborate)  are  higher.  However, 
BOH*  will  react  with  OH*  overcoming  0.55  eV  to  form  B(OH)2.  It  is 
possible  that  at  this  point  the  reaction  does  not  go  forward  due  to 
the  high  activation  barrier  between  B(OH)2  and  B(OH)3  and  the 
final  product  will  be  B(OH)2.  If  the  B-OH  activation  barrier  from 
B(OH)2  was  lower,  a  final  product  B(OH)3  might  react  in  the  solu¬ 
tion  phase  to  produce  the  hydrated  borate  anion  (B02)  [16]: 

B(OH)3(aq)  +  OH(-aq)  ~B(OH)4(aq)  (16) 

Then,  the  reaction  mechanism  would  be  finished  with  the 
release  of  8e~.  However,  further  investigations  are  required  to 
analyze  the  effect  of  the  solvating  conditions,  which  might  decrease 
the  activation  barriers  in  the  steps  involving  OH-.  Experimental 
techniques,  such  as  Raman  spectroscopy  would  also  help  to 
demonstrate  if  the  intermediate  B(OH)2,  is  especially  stable 
adsorbed  on  the  electrode  surface. 


Fig.  11.  Reaction  free  energy  diagram  for  borohydride  oxidation  over  Pd2Ir2(lll)  at  -0.5  V  vs.  SHE,  [BH4]  =  0.02  mol  dm-3,  [OH-]  =  2  mol  dm  3,  T  =  298  K.  The  red  line  represents 
the  most  favorable  mechanism  of  borohydride  oxidation  according  to  the  relative  energies  and  activation  barriers  calculated  from  DFT.  The  broken  and  dashed  lines  are  other 
possible  mechanisms  of  borohydride  oxidation.  Numbers  on  the  plot  label  the  activation  barriers  in  eV.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  this  article.) 
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Comparing  the  two  diagrams  shown  in  Figs.  10  and  11,  the 
oxidation  kinetics  are  more  favorable  on  Pd2l2  than  on  Pd2lir,  as  the 
energy  barriers  for  both  B-H  dissociation  and  B-OH  formation  are 
generally  lower.  The  Pd2lr2  alloy  catalyst  is  slightly  more  effective 
towards  the  borohydride  oxidation  reaction.  The  energy  barriers 
from  BOH"  to  B(OH)2  is  much  lower  on  Pd2lr2  than  on  Pd2lri,  which 
means  that  at  a  larger  concentration  of  Ir  in  the  alloy,  more  com¬ 
plete  oxidation  may  be  realized.  Other  than  that,  the  reaction 
mechanism  on  both  Pd-Ir  surfaces  follow  the  same  path,  for  which 
molecular  structures  are  summarized  in  Fig.  12. 

5.3.  Competitive  borohydride  oxidation  vs.  hydrolysis  at  the  Pd2\r\ 
(111)  surface 

During  the  borohydride  oxidation  on  Pd-Ir  surfaces,  hydrogen 
atoms  from  the  dehydrogenation  processes  are  adsorbed  on  the 
electrode  surface.  Those  hydrogen  atoms  might  either  combine 
with  another  adsorbed  hydrogen  atom,  forming  hydrogen  gas 
(Reaction  (17)),  or  oxidize  with  a  hydroxyl  ion  to  generate  water 
and  release  one  electron  (Reaction  (18)). 


H*  — >  1  /2H2  gas  +  * 

(17) 

H*  +  OHaq-H2Oaq  +  *  +  e- 

(18) 

The  free  energy  of  hydrogen  oxidation  and  hydrolysis  reactions 
can  be  calculated  by  using  Equations  (19)  and  (20),  respectively,  at 
1  atm  and  298  K.  A  borohydride  concentration  of  0.02  mol  dm'3 
and  3  mol  dm-3  NaOH  was  considered  to  calculate  the  free  energy 
of  H*  adsorption. 

AGEvol  =  0-5GH2gas  +  Qurf  ~  GH*  0^) 

AC^xid  =  Ch20  +  GSUrf  -  ^H*  -  COH“  -  eE  (20) 

For  an  accurate  consideration  of  the  competition  between 
hydrogen  evolution  and  hydrogen  oxidation,  the  hydrogen 
coverage  should  be  considered,  as  it  plays  an  important  role  in  the 
hydrogen  generation  rate.  The  hydrogen  coverage  is  dictated  by  the 
amount  of  surface  hydrogen  generated  by  borohydride 
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Fig.  12.  Configurations  of  the  preferred  borohydride  oxidation  path  on  Pd2Ir2(  111 )  surface.  The  dark  blue  atoms  correspond  to  Ir,  the  light  blue  big  atoms  correspond  to  Pd,  the  pink 
atoms  are  B,  the  red  are  O  atoms  and  the  white  are  H  atoms.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 
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dehydrogenation  and  removed  from  H2  desorption  or  H*  oxidation. 
For  a  given  voltage,  higher  H*  coverage  will  be  generated  for  low 
alkaline  solutions  (low  pH)  or  in  solutions  at  higher  temperatures. 
For  a  low  hydrogen  coverage  model,  one  hydrogen  atom  was 
optimized  on  the  Pd-Ir  surfaces,  giving  a  hydrogen  coverage  of  1  /9 
monolayer  (ML)  atom  for  the  3x3  Pd2-Iri  and  1/16  ML  atom  on 
the  4x4  Pd2-Ir2  surface.  For  comparison,  the  hydrogen  adsorption 
was  also  studied  on  a  Pd(lll)  surface.  For  that  purpose,  a  4-layer 
Pd(lll )  surface  was  optimized  in  a  3  x  3  cell.  The  lowest  energy  for 
the  optimized  hydrogen  atom,  among  all  the  possible  configura¬ 
tions  considered,  was  that  when  the  adsorbed  hydrogen  was  on  the 
fee  configuration  at  the  Pd(lll)  surface,  on  the  Pd-Ir  bridge  for 
Pd2lri( 111 )  and  on  the  top  of  the  Ir  atom  for  Pd2lr2(lll). 

The  comparison  of  the  free  energy  of  the  hydrogen  generation 
reaction  (Reaction  (18))  and  hydrogen  oxidation  reaction  (Reaction 
(17))  vs.  potential  for  low  hydrogen  coverage  on  Pd(lll),  Pd2-Iri 
and  Pd2-Ir2  surfaces  is  shown  in  Fig.  13.  The  free  energy  of 
hydrogen  evolution  did  not  change  with  the  concentration  of  Ir  in 
the  alloy  and  has  a  potential  independent  positive  value  of  0.43  eV. 
At  low  hydrogen  coverage,  hydrogen  evolution  is  not  favorable 
(AG  >  0)  for  either  Pd2-Iri  or  Pd2-Ir2  surfaces,  at  any  studied  po¬ 
tential.  The  free  energy  of  oxidation  was  positive  at  potentials 
from  -1  to  -0.4  V  vs.  Hg/HgO  (-0.86  to  -0.26  V  vs.  SHE),  sug¬ 
gesting  that  the  reaction  will  not  spontaneously  occur  if  the  H* 
coverage  is  low  at  those  potentials.  At  more  positive  potentials 
than  -0.4  V  vs.  Hg/HgO  (-0.26  V  vs.  SHE),  the  H*  oxidation  is 
favorable  (AG  <  0).  At  the  Pd(lll)  surface,  both  the  free  energy  of 
hydrogen  evolution  and  oxidation  were  more  positive  than  on  the 
Pd-Ir  surfaces.  Therefore,  at  low  H*  coverage  and  at  any  studied 
potential,  the  hydrogen  oxidation  and  evolution  are  more  favorable 
on  the  bimetallic  surfaces  than  on  pure  Pd(lll),  due  to  weaker  H 
binding  to  the  surface.  The  binding  energy  (BE)  of  the  H  to  the 
Pd(lll),  Pd2lri( 111 )  and  Pd2lr2(lll)  surfaces,  which  can  be  calcu¬ 
lated  using  Equation  (21)  [14  ,  are  -0.61  eV,  -0.43  eV  and  -0.41  eV, 
respectively. 

BE  —  £dft,h*  —  E*  —  Eh  (21 ) 

Weaker  H  binding  on  the  bimetallic  surfaces  might  lead  to  lower 
H  coverage,  and  therefore  a  faster  oxidation  reaction  rate  relative  to 
evolution,  as  the  oxidation  (rate  a  Oh)  and  evolution  rates  (rate  a 
0h)  directly  depend  on  the  hydrogen  coverage. 


Fig.  13.  Reaction  free  energies  for  H*  oxidation  (straight  line)  and  evolution  (dash  line) 
of  H2  gas  as  a  function  of  potential  for  low  hydrogen  coverage  1/9  ML  IT  on  Pd2lri(lll) 
surface  (black  lines),  1/16  ML  H*  on  Pd2Ir2( 111 )  surface  (red  lines)  and  1/9  ML  H*  on 
Pd(lll)  surface  (blue  lines).  The  Pd2Iri( 111 )  and  the  Pd2Ir2( 111 )  lines  overlay  each 
other.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 


At  higher  hydrogen  coverage,  the  free  energy  of  evolution  and 
oxidation  are  expected  to  become  more  exergonic,  as  increasing  the 
number  of  hydrogen  atoms  on  the  surface  destabilizes  H*.  At  high 
hydrogen  coverage,  9  hydrogen  atoms  were  optimized  on  the 
3x3  cell,  to  give  hydrogen  coverage  of  1  ML  atom  for  the  Pd2-Iri 
surface.  For  the  Pd2-Ir2  surface,  16  hydrogen  atoms  were  optimized 
on  the  4x4  cell  giving  1  ML  coverage.  The  large  number  of  possible 
configurations  of  the  9  hydrogen  atoms  on  the  surface  decreases 
the  likelihood  of  locating  a  global  optimum  structure.  Moreover, 
the  data  obtained  suffers  from  model  imperfections,  as  the  bime¬ 
tallic  structure  disrupts  the  symmetry  and  can  lead  to  effects  in  H 
binding  that  are  not  practical  in  a  typical  experimental  alloy  elec¬ 
trode,  which  is  probably  less  ordered.  The  difference  in  cell  size 
between  Pd2lri  and  Pd2lr2  also  allows  for  differences  in  the  surface 
reorganization  allowed.  However,  these  calculations  can  provide  an 
orientation  of  the  hydrogen  behavior  at  high  coverage  on  the  Pd-Ir 
and  Pd  surfaces. 

At  the  pure  Pd  electrode,  the  optimized  structure  had  the  9H 
atoms  in  fee  sites.  In  contrary,  on  the  Pd2lri  surface,  the  9H  atoms 
changed  location  during  the  optimization  process  from  the  sym¬ 
metric  hollow  sites  to  interact  with  the  Ir  atoms.  At  the  Pd2lr2 
surface,  the  16  hydrogen  atoms  were  initially  located  on  Pd  and  Ir 
atop  sites,  but  the  hydrogen  atoms  tend  to  move  towards  the  Ir 
instead  of  to  the  Pd  atoms.  The  Pd2lr2  surface  to  offers  enough  Ir 
atoms  such  that  all  H  can  move  from  their  symmetric  hollow  sites 
towards  Ir  atoms,  making  stripes  of  H  atoms.  These  Ir  stripes,  which 
form  a  regular  pattern,  might  provide  a  greater  stabilization  to  the 
H  high  coverage  and  the  strong  H  binding  could  hurt  the  selectivity 
of  the  catalyst  towards  the  oxidation.  However,  these  calculations 
consider  an  ideal  material  configuration  consisting  of  Pd  layers 
followed  by  Ir  layers  that  are  difficult  to  achieve  experimentally. 

The  results  showed  more  negative  and  thus  more  favorable  free 
energies  of  hydrogen  oxidation  and  hydrogen  evolution  than  those 
obtained  at  low  hydrogen  coverage,  as  was  expected.  At  high 
hydrogen  coverage  on  the  Pd2lri  surface,  the  free  energy  of  the 
hydrogen  evolution  was  favorable  (approximately  -0.73  eV).  That 
value  was  0.7  eV  more  negative  on  Pd2lr2  suggesting  that  the 
presence  of  Ir  makes  more  favorable  the  hydrogen  generation.  The 
crossing  point  at  which  oxidation  competes  with  evolution  does 
not  vary  for  the  studied  materials  (Pd,  Pd2lri  and  Pd2lr2).  The  two 
lines  in  Fig.  13,  corresponding  to  free  energy  of  hydrogen  evolution 
and  oxidation,  intercept  at  approximately  -1  V  vs.  Hg/HgO  (-0.86  V 
vs.  SHE),  and  the  free  energy  of  oxidation  is  lower  than  the  free 
energy  of  hydrogen  evolution  at  more  positive  potentials.  That 
suggests  that  at  more  positive  potentials  than  -1  V  vs.  Hg/HgO 
(-0.86  V  vs.  SHE),  the  hydrogen  atom  adsorbed  on  the  surface 
would  oxidize  rather  than  evolve  as  hydrogen  bubbles.  At  more 
negatives  potentials  than  -1  V  vs.  Hg/HgO  (-0.86  V  vs.  SHE),  the 
hydrogen  evolution  is  more  favorable  than  the  oxidation  and  thus, 
hydrogen  bubbles  would  evolve.  This  agrees  with  the  experimental 
results  showed  in  Fig.  2,  as  hydrogen  generation  was  observed  at 
the  open  circuit  potential  (-0.99  V  vs.  Hg/HgO  or  -0.85  V  vs.  SHE) 
and  at  more  positive  potentials  the  hydrogen  generation  substan¬ 
tially  decreased.  It  can  be  concluded  that  the  presence  of  Ir  in  the 
bimetallic  alloy  decreases  the  free  energy  of  oxidation,  making  the 
Pd2lr2  surface  more  favorable  for  the  borohydride  oxidation  at  all 
potentials,  compared  to  the  Pd2lri. 

6.  Summary  and  conclusions 

The  catalytic  activity  of  the  Pd-Ir  alloy  for  the  borohydride 
oxidation  and  its  hydrolysis  was  experimentally  and  computa¬ 
tionally  analyzed.  The  Pd-Ir  coated  microfibrous  carbon  supported 
on  a  titanium  plate  electrode  showed  catalytic  activity  towards  the 
borohydride  oxidation,  obtaining  current  densities  from  100  to 
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200  mA  cm-2,  when  a  solution  containing  0.5  mol  dm-3  NaBH4  in 
3  mol  dm”3  NaOH  and  a  3  cm2  electrode  geometric  surface  area, 
were  used.  At  the  same  time,  the  hydrogen  generation  rate  was 
minimized  to  less  than  0.1  dm3  min-1.  For  that  reason,  Pd-Ir  seems 
to  be  a  promising  material  for  the  borohydride  oxidation.  However, 
there  are  some  drawbacks  that  need  to  be  addressed  before  the 
Pd-Ir  alloy  can  be  used  in  a  DBFC,  such  as  the  high  cost  and  the 
spontaneous  hydrogen  generation  at  open  circuit  potential 
(-0.99  V  vs.  Hg/HgO,  -0.85  V  vs.  SHE).  Further  experiments  need  to 
be  carried  out  in  a  complete  cell,  in  order  to  test  the  catalyst  and 
estimate  the  importance  of  the  inconvenience  of  the  hydrogen 
generation  at  the  OCP. 

DFT  studies  were  carried  out  at  the  Pd2lri  and  Pd2lr2  surfaces  in 
order  to  analyze  the  mechanism  of  the  borohydride  oxidation  re¬ 
action  and  evaluate  the  effect  of  the  Ir  concentration.  The  results 
showed  that  in  the  presence  of  Ir,  the  Gibbs  free  energy  of  BH4 
adsorption  is  lower  than  that  on  pure  Pd,  at  all  potentials,  which 
suggests  that  the  oxidation  reaction  starts  at  more  negative  po¬ 
tentials  on  the  alloy  than  on  pure  Pd.  DFT  studies  also  showed  that 
at  high  hydrogen  coverage  and  -1  V  vs.  Hg/HgO  (-0.86  V  vs.  SHE), 
there  is  hydrogen  generation  but  at  more  positive  potentials  the 
oxidation  of  borohydride  is  more  favorable  than  the  hydrolysis. 
These  results  agree  with  experimental  observation,  the  reaction 
mechanism  on  both  Pd-Ir  surfaces  followed  the  same  path, 
showing  a  favorable  8  electron  transfer.  However,  the  activation 
energies  became  too  high  in  the  last  elementary  steps,  which 
involve  B-OH  binding.  The  activation  barriers  were  generally  lower 
on  Pd2lr2  than  on  Pd2lri,  with  the  reaction  mechanism  more 
favorable  on  the  former.  From  BOH  at  Pc^In  and  BOHOH  at  Pd2Ir2, 
the  activation  barriers  became  too  high  and  it  is  possible  that  the 
reaction  does  not  continue  any  further.  In  that  case,  the  number  of 
electrons  released  would  be  5  and  6,  at  Pd2lri  and  Pd2lr2,  respec¬ 
tively.  Those  intermediate  species,  BOH  at  Pd2lri  and  BOHOH  at 
Pd2lr2,  might  cover  the  catalyst  surfaces,  being  a  poison  and 
reducing  the  catalytic  activity  of  the  electrodes  towards  the  boro¬ 
hydride  oxidation. 

Although  the  complete  oxidation  reaction  seems  to  be  favorable 
according  to  the  relative  energy  diagram,  it  is  possible  that  the  DFT 
approach  used  for  the  calculation  of  the  activation  barriers  of  the 
B-OH  forming  reactions  overestimates  energy  barriers,  due  to  the 
lack  of  consideration  of  solvation.  Computational  models  that 
include  the  role  of  water,  by  adding  two  H2O  molecules  next  to  the 
BH4  molecule,  can  be  utilized  to  decrease  the  activation  barriers 
corresponding  to  B-OH  bond  formation  and  obtain  the  final 
oxidation  product,  metaborate,  and  release  8e-. 


References 

[1]  C.  Ponce  de  Leon,  F.C.  Walsh,  D.  Pletcher,  D.J.  Browning,  J.B.  Lakeman,  J.  Power 
Sources  155  (2006)  172-181. 

[2]  I.  Merino-Jimenez,  C.  Ponce  de  Leon,  A.  Shah,  F.C.  Walsh,  J.  Power  Sources  219 
(2012)  339-357. 

[3]  B.H.  Liu,  Z.P.  Li,  S.  Suda,  Electrochim.  Acta  49  (2004)  3097-3105. 

[4]  M.V.  Mirkin,  H.  Yang,  A.J.  Bard,  J.  Electrochem.  Soc.  139  (1992)  2212—2217. 

[5]  S.C.  Amendola,  P.  Onnerud,  M.T.  Kelly,  P.J.  Petillo,  S.L.  Sharp-Goldman, 
M.  Binder,  J.  Power  Sources  84  (1999)  130—133. 

[6]  B.M.  Concha,  M.  Chatenet,  E.A.  Ticianelli,  F.H.  Lima,  R.D.  Lima,  216  ECS  Meeting 
Abstracts  11th  International  Symposium  on  Solid  Oxide  Fuel  Cells,  vol.  902, 
Vienna,  Austria,  2009,  p.  368. 

[7]  M.  Chatenet,  F.H.B.  Lima,  E.A.  Ticianelli,  J.  Electrochem.  Soc.  157  (2010) 
B697-B704. 

[8]  E.  Gyenge,  M.  Atwan,  D.  Northwood,  J.  Electrochem.  Soc.  153  (2006) 
A150-A158. 

[9]  J.  Ma,  N.A.  Choudhury,  Y.  Sahai,  Renewable  Sustainable  Energy  Rev.  14  (2010) 
183-199. 

[10]  D.A.  Finkelstein,  N.D.  Mota,  J.L.  Cohen,  H.c.D.  Abruria,  J.  Phys.  Chem.  C  113 
(2009)  19700-19712. 

[11]  G.H.  Miley,  N.  Luo,  J.  Mather,  R.  Burton,  G.  Hawkins,  L.  Gu,  E.  Byrd,  R.  Gimlin, 
P.J.  Shrestha,  G.  Benavides,  J.  Laystrom,  D.  Carroll,  J.  Power  Sources  165  (2007) 
509-516. 

[12]  J.Q,  Yang,  B.H.  Liu,  S.  Wu,  J.  Power  Sources  194  (2009)  824-829. 

[13]  V.  Kiran,  T.  Ravikumar,  N.T.  Kalyanasundaram,  S.  Krishnakumar,  A.K.  Shukla, 
S.  Sampath,  J.  Electrochem.  Soc.  157  (2010)  B1201— B1208. 

[14]  G.  Rostamikia,  M.J.  Janik,  Electrochim.  Acta  55  (2010)  1175—1183. 

[15]  G.  Rostamikia,  M.J.  Janik,  J.  Electrochem.  Soc.  156  (2009)  B86— B92. 

[16]  G.  Rostamikia,  A.J.  Mendoza,  M.A.  Hickner,  M.J.  Janik,  J.  Power  Sources  196 
(2011)  9228-9237. 

[17]  J.K.  Norskov,  J.  Rossmeisl,  A.  Logadottir,  L.  Lindqvist,  J.R.  Kitchin,  T.  Bligaard, 
H.  Jonsson,  J.  Phys.  Chem.  B  108  (2004)  17886-17892. 

[18]  X.  Nie,  M.R.  Esopi,  M.J.  Janik,  A.  Asthagiri,  Angew.  Chem.  Int.  Ed.  52  (2013) 
2459-2462. 

[19]  C.J.  Patrissi,  R.R.  Bessette,  Y.K.  Kim,  C.R.  Schumacher,  J.  Electrochem.  Soc.  155 
(2008)  B558-B562. 

[20]  C.  Celik,  F.G.  Boyaci  San,  H.I.  Sarac,  J.  Power  Sources  185  (2008)  197—201. 

[21]  B.H.  Liu,  Z.P.  Li,  S.  Suda,  J.  Electrochem.  Soc.  150  (2003)  A398-A402. 

[22]  K.  Wang,  J.  Lu,  L.  Zhuang,  Catal.  Today  170  (2011)  99-109. 

[23]  A.E.  Bolzan,  J.  Electroanal.  Chem.  380  (1995)  127—138. 

[24]  G.  Denuault,  C.  Milhano,  D.  Pletcher,  Phys.  Chem.  Chem.  Phys.  7  (2005) 
3545-3551. 

[25]  Z.X.  Liang,  T.S.  Zhao,  J.B.  Xu,  L.D.  Zhu,  Electrochim.  Acta  54  (2009)  2203—2208. 

[26]  A.B.  Anderson,  D.B.  Kang,  J.  Phys.  Chem.  A  102  (1998)  5993—5996. 

[27]  G.  Kresse,  J.  Furthmiiller,  Comput.  Mater.  Sci.  6  (1996)  15—50. 

[28]  G.  Kresse,  J.  Furthmiiller,  Phys.  Rev.  B  54  (1996)  11169. 

[29]  G.  Kresse,  J.  Hafner,  Phys.  Rev.  B  47  (1993)  558. 

[30]  A.U.  Nilekar,  M.  Mavrikakis,  Surf.  Sci.  602  (2008)  L89— L94. 

[31]  Y.  Xu,  A.V.  Ruban,  M.  Mavrikakis,  J.  Am.  Chem.  Soc.  126  (2004)  4717-4725. 

[32]  J.  Greeley,  T.F.  Jaramillo,  J.  Bonde,  I.  Chorkendorff,  J.K.  Norskov,  Nature  5 
(2006)  909-913. 

[33]  H.J.  Monkhorst,  J.D.  Pack,  Phys.  Rev.  B  13  (1976)  5188-5192. 

[34]  G.  Helkeman,  B.P.  Uberuaga,  H.  Jonsson,  J.  Chem.  Phys.  113  (2000)  9901. 

[35]  A.J.  Bard,  L.R.  Faulkner,  Electrochemical  Methods  Fundamentals  and  Applica¬ 
tions,  second  ed.,  John  Wiley  and  Sons,  Inc.,  2001. 

[36]  D.  Tsiplakides,  C.G.  Vayenas,  Solid  State  Ionics  152—153  (2002)  625—639. 

[37]  G.  Rostamikia,  R.J.  Patel,  I.  Merino-Jimenez,  M.  Hickner,  M.J.  Janik,  in 
preparation. 


